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ABSTRACT: In the context of the decoherent histories approach to the quan- 
tum mechanics of closed systems, Gell-Mann and Hartle have argued that the vari- 
ables typically characterizing the quasiclassical domain of a large complex system 
are the integrals over small volumes of locally conserved densities - hydrodynamic 
variables. The aim of this paper is to exhibit some simple models in which ap- 
proximate decoherence arises as a result of local conservation. We derive a formula 
which shows the explicit connection between local conservation and approximate 
decoherence. We then consider a class of models consisting of a large number 
of weakly interacting components, in which the projections onto local densities 
may be decomposed into projections onto one of two alternatives of the individual 
components. The main example we consider is a one-dimensional chain of locally 
coupled spins, and the projections are onto the total spin in a subsection of the 
chain. We compute the decoherence functional for histories of local densities, in 
the limit when the number of components is very large. We find that decoherence 
requires two things: the smearing volumes must be sufficiently large to ensure ap- 
proximate conservation, and the local densities must be partitioned into sufficiently 
large ranges to ensure protection against quantum fluctuations. 
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1. INTRODUCTION 



One of the primary aims of quantum cosmology is to understand how the uni- 
verse is classical to a very high degree of precision, given the hypothesis that it 
is described at the most fundamental level by quantum theory [1,2,3]. Mathe- 
matically, this aim translates into the question of why, on sufficiently large scales, 
quantum theory admits an emergent description of the universe involving only a 
small number of dynamical variables obeying an approximately closed set of deter- 
ministic evolution equations. 

Whilst many approaches to this question concentrate on the positions and 
momenta of point particles, these dynamical variables are special cases of a more 
general description of the physical world in terms of local densities: number density, 
momentum density, energy density, charge density, etc. The question of emergent 
classicality then consists first, of understanding why these variables enjoy such a 
distinguished role, and second, of deriving the familiar hydrodynamic equations for 
these densities. 

The decoherent histories approach to quantum theory is a recently devel- 
oped formulation of quantum theory which is particularly suited to this problem 
[4,5,6,7,8,9]. In brief, the approach permits predictions to be made in genuinely 
closed systems, such as the entire universe, without relying on notions of measure- 
ment or on an external classical domain. The goal of the approach is to assign 
probabilities to the histories of a closed quantum system. It is this feature that 
makes it particularly useful for discussing emergent classicality. This is because to 
assert in quantum theory that a certain variable approximately satisfies a deter- 
ministic evolution equation involves computing the probability for a time-ordered 
sequence of values of that variable, ie., for a history. It is then necessary to show 
that the probability for that history is strongly peaked about the sequence of values 
corresponding to the evolution equation. 

This paper constitutes the beginnings of a general application of the decoherent 
histories approach to the question of deriving hydrodynamic equations. 

We begin with a very brief review of the decoherent histories approach [4,5,6,7,8,9,10] , 
followed by a general discussion of the issue of emergent classicality for hydrody- 
namic variables. A more detailed statement of the particular models considered in 
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this paper may be found at the end of this section. 

The models and technical results presented here are admittedly simple, perhaps 
to the extent of being overshadowed by the very broad discussion of emergent 
classicality given below. This paper is, however, a small first step in a possibly 
extensive programme. We have therefore chosen to give a general sketch of that 
programme, to indicate the direction of future research and to provide the general 
context in which the possible significance of our modest results may be understood. 

1(A). The Decoherent Histories 
Approach to Quantum Theory 

A quantum-mechanical history is a sequence of propositions at a succession of 
times. Propositions at a fixed moment of time are represented by sets of projections 
operators, {Pa}- They are exhaustive and exclusive, 

^Pa = l, PaPp = S^pPa (1.1) 
a 

A projector is said to be fine-grained if it is of the form |q;)(q;|, where {\a)} are 
a complete set of states. Otherwise it is said to be coarse-grained. A single 
quantum-mechanical history is characterized by a string of time-dependent projec- 
tions, P^^{ti), ■ ■ ■ ,P^^{tn), together with an initial state p. The time-dependent 
projections are related to the time-independent ones by 

Paki^k) = e'^^^^-^^^P^^e-'^^*''-^^^ (1.2) 

where H is the Hamiltonian. The candidate probability for such histories is 

p{ai,a2, ■■■an) = Tr [P^Jin) ■ • •<(tl)p<(tl) • ■ ■ PSJtn)) (1-3) 

To be a true probability, Eq.(1.3) must satisfy the probability sum rules. That 
is, it must be such that the probability for each coarser-grained history is the 
sum of the probabilities for the constituent finer-grained histories. Coarser-grained 
histories may be constructed, for example, by summing the projections at each 
moment of time: 

Pa=J2P^ (1-4) 
aGa 
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The probability sum rules are generally not satisfied due to quantum interference. 

Sets of histories which do not sufi^er interference, and hence for which the sum 
rules are satisfed may be found using the decoherence functional, 

D{a,^) = Tv [P^^tn) ■ ■ ■ Pi,{ti) pP\{ti) ■ ■ ■ P^>^{tn)) (1.5) 

Here a denotes the string ai,a2, ■ ■ - (Xn- Intuitively, the decoherence functional 
measures the amount of interference between pairs of histories. It may be shown 
that the probability sum rules are satisfied for all coarse-grainings if and only if 

ReD{a,a) = (1.6) 

for all distinct pairs of histories a, a' [7]. Such sets of histories are said to be 

consistent (or weakly decoherent). 

The consistency condition (1.6) is usually (but not always) satisfied only for 
coarse-grained histories. When sets of histories satisfy the consistency condition 
(1.6) as a result of coarse-graining, they often satisfy, in addition, the stronger 
condition of decoherence, 

£>(«,«') =0, for a^a' (1.7) 

(sometimes classified as medium decoherence, with yet more stringent criteria for 
strong decoherence [6]). Physically, decoherence is intimately related to the exis- 
tence of records about the system somewhere in the universe [5] . 

In most cases of interest decoherence is only approximate, so measures of ap- 
proximate decoherence are required. We say that the degree of decoherence is of 
order e if the probability sum rules are violated only up to terms of order e times 
the probabilities themselves [10]. It may be shown that, under certain conditions 
on the distribution of off-diagonal terms in the decoherence functional, decoherence 
to order e will be satisfied if 

\D{a,g!)\'^ < €^ D{a,a)D{a',a'). (1.8) 

1(B). Emergent Quasiclassicality 

We would like to use the decoherent histories approach to demonstrate the 
emergence of an approximately classical world from an underlying quantum one, 
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together with the quantum fluctuations about it described by the standard Copen- 
hagen quantum mechanics of measured subsystems. Such a state of affairs is re- 
ferred to as a quasiclassical domain [4,5,6]. In more technical terms, a quasiclassical 
domain consists of a decoherent set of histories, characterized largely by the same 
types of variables at different times, and whose probabilities are peaked about 
deterministic evolution equations for the variables characterizing the histories. 

The histories should, moreover, be maximally refined with respect to a specified 
degree of approximate decohcrcnce. That is, one specifies a decoherence factor 
e in the approximate decohcrcnce condition discussed above. This should, for 
example, be chosen so that the probabilities are defined to a precision far beyond 
any conceivable test. Then, the histories should be fine-grained {e.g., by reducing 
the widths of the projections) to the point that further fine-graining would lead to 
violation of the specified degree of approximate decoherence. The resulting set of 
histories are then called maximally refined. The reason for maximally refining the 
histories is to reduce as much as possible any apparent subjective element in the 
choice of coarse-graining. 

Given the Hamiltonian and initial state of the system, one's task is to compute 
the decoherence functional for various different choices of histories, and see which 

ones lead to quasiclassical behaviour. We expect this to be a formidable task, and 
there is no reason to believe that there will be a unique answer. Generally, one 
might expect that there will be a hierarchy of variables, similar to the BBGKY 
hierarchy [11]. 

Many previous discussions [5,10,12,13,14] of emergent classicality concern sys- 
tems in which there is a natural separation of the total closed sytem into "system" 
and "environment" , and this separation is the source of the coarse-graining required 
for decoherence. t A generic closed system, however, will usually not have such a 
separation, and it is one of the strengths of the decoherent histories approach that 
it does not rely on the existence of one. Gertain variables will, however, be distin- 
guished by the existence of conservation laws for total energy, momentum, charge, 
particle number, etc. Associated with such conservation laws are local conservation 

t More precisely, by separation into system and environment, we mean that the total Hilbert 
space for the closed system may be written as a tensor product of the system and envi- 
ronment Hilbert spaces. Some authors appear to use this expression to mean something 
more general. 
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laws of the form 

|+V.J = (1.9) 

The candidate quasiclassical variables are then 

Qv= [ d^x p{x) (1.10) 

Jv 

These are the local densities discussed above. If the volume V becomes infinite, 
Qy will be an exactly conserved quantity. In quantum mechanics it will commute 
with the Hamiltonian, and, as is easily seen, histories of Qy's will decohere exactly 

[15] 

If the volume is finite but large compared to the microscopic scale, Qy will be 
slowly varying compared to other dynamical variables. This is because the local 
conservation law (1.9) permits Qy to change only by redistribution, which is limited 
by the rate at which the locally conserved quantity can fiow out of the volume. 
Because these quantities are slowly varying, histories of them should approximately 
decohere. 

More precisely, when the volume V is infinite in size, the decoherence functional 
for histories of the variables (1.10) will be exactly diagonal. As V shrinks from infi- 
nite size, the decoherence functional will develop off'-diagonal terms. However, one 
would not expect the off-diagonal terms to grow appreciably until V approaches 
the length scales characteristic of the system at hand. These length scales can de- 
pend only on the initial state and on the Hamiltonian. The length scale associated 
with the initial state could in principle take any value, hence some restriction on 
the initial state will be necessary. The length scales associated with Hamiltonians 
of the type we are typically interested in will generally be very small. The various 
physically relavent scales, are, for example, the mean free path between molecular 
collisions, or the length scale of intermolecular forces. Therefore, for some class of 
initial states (to be determined), one would expect the off-diagonal terms of the 
decoherence functional to remain small until V shrinks down to the microscopic 
scale. 

Hence, for a suitable class of initial states, we expect the variables (1.10) to be 
approximately decoherent as a consequence of their association with conservation 
laws, which are in turn connected with the global symmetries of the system. 

Given decoherence, we may then examine the probabilities for histories of hy- 
drodynamic variables, and ask whether they are peaked about the expected hydro- 
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dynamic equations. Derivations of the hydrodynamics equations from an underly- 
ing microscopic quantum theory have certainly been carried out before (see Ref.[16] 
for example). These derivations have shown that the expectation values of the lo- 
cal densities Qy, in a local equilibrium state, evolve in approximate accordance 
with the hydrodynamic equations. The derivation contemplated here, however, is 
considerably more general (although undoubtedly related in some way). In the 
decoherent histories approach we would like to show that the probabilities for his- 
tories of imprecisely specified values of the local densities, for a more general class 
of initial states, are peaked about the hydrodynamic equations. This paper is only 
a first step in that direction; but it is important to keep the ultimate goal in mind. 

The connection between the earlier derivations of the hydrodynamic equations 
and the histories derivation contemplated here is analagous to the connection be- 
tween the density operator and the decoherent histories approaches to emergent 
classicality of quantum Brownian motion models. In the density operator approach, 
it was shown that an initial Gaussian wave packet followed an approximately clas- 
sical path, with dissipation, and with fluctuations due to classical and quantum 
noise [13,12,17]. In the decoherent histories approach, it was shown that for a 
wide variety of initial states, the probabilities for histories of position samplings is 
strongly peaked about classical evolution, with dissipation, with a width of peaking 
depending on classical and quantum noise [5,10]. 

1(C). This Paper 

The programme sketched above is clearly a very extensive one. Although clear 
in principle how to proceed, it is very difficult to carry out in practice. In principle, 
since the operators are known, one may compute their spectrum, from which the 
projectors onto ranges of the spectrum may be derived. The time evolution of 
those projectors may be deduced and inserted in the decoherence functional. In 
practice, this is very difficult and it seems likely that new mathematical techniques 
will be needed. 

It is unlikely, for example, that the path integral techniques so successful in 
the study of quantum Brownian motion models will be eff^ective here [5,10]. In the 
quantum Brownian motion models, the interesting variables are the position or mo- 
mentum of a distinguished particle, and projections onto these variables are easily 
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implemented as restrictions on the paths in a sum over paths. The hydrodyanmic 
variables considered here, however, are non-trivial functions of positions and mo- 
menta, and projections onto ranges of their spectra cannot in general be expressed 
in terms of restrictions in a sum over paths in phase space. 

The quantum Brownian motion models also made heavy use of the influence 
functional method [18,19]. Again this is inapplicable, because it relies on being 
able to explicitly integrate out the "environment" {i.e., the variables ignored in 
the coarse-graining procedure), which cannot be done here. 

To make a start on the general problem what is required is some very simple 
models which retain the essential features of approximate decoherence through 
approximate conservation, yet are simple enough to be solvable. In this paper we 
will present some models which are of this type. 

In Section II, we derive a formula expressing the general connection between 
approximate conservation and approximate decoherence. We will not use this for- 
mula explicitly in this paper, although we present it to flesh out the heuristic 
arguments given above, and because it may be of practical use elsewhere. 

In Section III, we consider a system consisting of a large number N of weakly 
interacting components in which the individual components are described by just 
two alternatives. An example is a large number of non-interacting particles in a box 
divided into two sections, and the two alternatives for each individual component 
are that the particle is in the left or right section of the box. A crude hydrodynamic 
variable is then the number of particles in, say, the right hand section. The pro- 
jection onto such a hydrodynamic variable can be given in terms of the projections 
onto the individual particles, and we derive a formula which gives this connection 
precisely. Another example is the spin system model we describe below. For this 
class of models, we compute the decoherence functional for histories characterized 
by projections onto densities of large collections of particles. We thus compute the 
degree of decoherence. 

The calculation of Section III is principally concerned with the combinatoric 
aspects of large collections of particles. The feature of decoherence through approx- 
imate conservation, with which we are ultimate concerned, is to be found in the 
decoherence functionals for the individual components. Hence we need to compute 
such a decoherence functional explicitly. This is carried out for a spin model in 

9 



Section IV. 

The model consists of a long chain of M locally coupled spins, and the Hamil- 
tonian conserves the total spin. We employ the simplest coarse-graining pos- 
sible, which is to divide the chain into two pieces, of lengths Mi and M2 (so 
Ml + M2 = M), and project onto the amount of spin into one section of the chain. 

The ground state of the system is the state in which all the spins point up. The 
first excited states, the so-called spin waves, consist of superpositions of particle-like 
states in which one of the spins is down. The higher excited states are more difficult 
to construct because they involve interactions between spin waves, but they may 
be approximated by suitable combinations of the first excited states in the limit 
that the coupling is weak. In this approximation the system is therefore not unlike 
a collection of weakly interacting particles. We thus show that the decoherence 
functional for projections onto the number of, say, down spins in one section may 
be approximated by a product of decoherence functionals for the "one particle" 
states, in which the projections ask in which side is the down spin. This means 
that, in this approximation, our spin model is a system of the type considered in 
Section III. 

We calculate the decoherence functional for the individual components numer- 
ically. It is approximately diagonal, but not surprisingly, the degree of decoherence 
is not very good, being no better than would be expected from the overlap of a 
random pair of states in a large Hilbert space. Its significant feature is the way 
the degree of decoherence changes as Mi and M2 are changed. We use this re- 
sult, together with the formula derived in Section III, to compute the degree of 
decoherence for projections onto ranges of spin density. 

We summarize and conclude in Section V. 
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2. APPROXIMATE CONSERVATION 
AND APPROXIMATE DECOHERENCE 



In this section we briefly outline the connection between approximate conser- 
vation and approximate decoherence. 

As in Section I, let Q denote a local density averaged over a volume V: 

Q = [ (fix p{^) (2.1) 
Jv 

By integrating the conservation equation for p, Eq.(1.9), over a time interval [0,t\ 
and over the spatial volume V, it is readily shown that 

Qt = Q-[*dt'f (fixn-j{yi,t') (2.2) 
Jo JdV 

where dV is the boundary of V. This relation will also hold in the quantum theory 
for suitably ordered operators. 

Now consider the decoherence functional for histories characterized by projec- 
tions onto these operators at two moments of time: 

D{ai, a2\a[, ^2) = Tr [Pa^e-'^* Pa.pPa'/^') (2.3) 

This may be written 

D(ai,a2|ai,«2) = Y.^m\Pa^e-'^^ Pc,\n) pnm{m\Pa^e-'^^ P^,^\nY (2.4) 
m,n 

where |n) denotes a complete set of states, and pnm — {n\p\m). The degree of 
decoherence will therefore depend on the size of the amplitudes, 

Amn{oci,a2) = {m\Pa^e-'^^ PaM) (2.5) 

Projections onto exactly conserved quantities will commute with H, and Amn{cti, ct2) 
will be exactly zero unless q;2 = ai, hence the decoherence functional will be ex- 
actly diagonal. Wc arc interested, however, in the case in which the projections 
commute only approximately with the Hamiltonian. 

Suppose the spectrum of the operator Q is continuous. Then projections onto 
precisely specified values of its eigenvalues may be constructed using the delta- 
function: 

Pa = 5{Q - a) (2.6) 
11 



(These are not fine-grained projections, since the eigenvalues will typically be highly 
degenerate). Projections onto imprecisely specified ranges of eigenvalues are readily 
constructed by summing over ranges of a, although we will not do this explicitly 
here. 

The amplitudes Amnictl^ ct2) iiow be written 

Amn{ai,a2) = {m\S{Q - e'^^ 5{Q - ai)\n) 

= {m\S{Q - a2) e-'^^ e'^^ e'^* S{Q - ai) \n) (2.7) 

where e is an arbitrary parameter. Now we use the fact that 

^icQt ^ ^iHt ^leQ ^-iHt (2.8) 

where Qt = e^^^Qe~^^^ . Hence we have 

Amn{ai,a2) = {m\S{Q - ^3) e"^^^ e"*^^ e*^^* S{Q - ai)\n) (2.9) 

Qt is also given by (2.2), so inserting (2.2), expanding to first order in e, using the 
fact that QS{Q — a) — ad{Q — ct), and rearranging, we get 

Amn{ai,a2) = r/ dt'f d^x {m\6{Q - a2) e'^^ n ■ 3{:)c,t') 6{Q - ai)\n) 

[a2 - ai) Jq Jqy 

(2.10) 

Eq.(2.10) is the main result of this section. Inserted in the decoherence functional 
Eq.(2.4), (2.10) yields the degree of decoherence as a function of the size of the 
smearing volume V. In particular, we see that decoherence becomes exact as the 
boundary dV of V goes to infinity, as it must, since the operators Q are then exactly 
conserved. Of course, the rate at which approximate decoherence approaches exact 
decoherence will depend on the states |n), |m), in accordance with our general 
expectations. 

In the case of exact conservation the amplitude Amn{ctl,ct2) is non-zero only 
for q;2 = cti, hence the probability of 0:2 given cti is 1 when 0:2 = cti, and zero oth- 
erwise. In the case of approximate conservation, however, one expects the situation 
to be a little more complicated. Given decoherence, the probability distribution 
for histories will not necessarily be peaked about 02 = cti. Rather, since we expect 
the probabilities for histories to indicate hydrodynamic equations, it will be peaked 
about values of ai, a2 consistent with some deterministic evolution equations. For 
this reason, it is perhaps better to speak of approximate determinism, rather than 
approximate conservation. 
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3. DECOHERENCE OF DENSITIES OF LARGE 
COLLECTIONS OF NON-INTERACTING PARTICLES 



3(A). The System 



Consider a system which consists of a very large number of particles or 
components whose interactions are so weak that they may be neglected. We are 
interested in the case in which each component is described by just two alternatives 
at each moment of time, which may be represented by projections P and P = 1 — P. 
An example is a collection of particles in a box divided into two sections, and 
the projections P and P then represent the propositions that the particle is in, 
respectively, the right-hand or left-hand section of the box. A version of this 
system involving spins will be described in the next section. We will refer to the 
alternatives represented by P and P as "yes" and "no" , respectively. 

Using these elementary projections onto the individual components of the sys- 
tem, projections onto densities of the whole system may be constructed. For a 
system of two particles, for example, the number of particles in the right-hand sec- 
tion of the box, may be 2, 1 or 0. These propositions are represented, respectively, 
by the projections. 



It is easy to see that these projections are mutually exclusive and exhaustive, as 
required. 

For large systems, the projections onto densities rapidly become quite cumber- 
some. However, the following trick turns out to be extremely useful. We have the 
identity. 



where 6nm is the Kronecker delta. The projection operator onto number density n 
in a system of N components is then given by 



P2 = P®P 




Pq = P®P 



Pl=P^P + P^P 



(3.2) 
(3.3) 




(3.4) 




dX e 



—iXn 



Fi(A) ®F2(A)---®F^(A) 



(3.5) 



TT 
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Here 

Ffc(A) = e*^p(^) + P(^) (3.6) 

where p(^) and p(^) denote the projections P, P operating on component number 
k. What happens in Eq.(3.5) is that in the tensor product over the F)^(A)'s, the 
projection onto number density n occurs with coefficient, e^"^", z.e., 

Fi(A) ® F2{\) ■■■® Fn{X) =P P • • • ® P 

+ e^'^(P(8)P<g)---(g)P + P(8)P(8)---(8)PH ) 

+ ••• 

N 

= J2 e^^" Pn (3.7) 
n=0 

The integration over A then picks out only the term with coefficent e*^"^. 

These projections are onto precise vahics of number density. Generally one is 
interested in imprecise values of number density, i.e., whether the number density 
lies in a specified range, and these projections are obtained by summing over n. 
We will consider this in more detail below. 

3(B). The Decoherence Functional 

We may now write down the decoherence functional for histories characterized 
by projections at two moments of time onto precise values of number density. It 
is, 

D{ni,n2\n[,n2) = TV [Pn^e''^^ Pn.pPn'/"* Pn,) 

(27r)* i-TT J-TT J-TT J-TT 

X Tr (Pi(A2) ® P2(A2) ■ ■ ■ ® Fn{X2) e''^^* Pi(Ai) ® P2(Ai) Fn{Xi) 
X pT f}{X[) ® pJ'(A'i) • • • ® Fi(A'i) e'"^' f}{X'2) ^ pJ(A'2) • • • ® pJ;r(A'2)) (3-8) 

Note that we can of course use the simplifying feature (Pfi,) = Pfi, the final 
projection, but it will become clear below why we have not done this. 

To make further progress we make two simplifications. First, we assume that 
the interaction between component particles is negligible, so the Hamiltonian Hx 
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for the total system has the form 

Ht = H®1---®1 + 1^H---®1 + --- (3.9) 

Second, we assume that the intial density operator p^p for the total system factors, 
and that the density operator for each component system is the same: 

pj' = P®p---®P (3.10) 
The decoherence functional may now be written 

[ZTTP J-TT J-TT J -IT 

N 

X n ^ (^^(^2) e-^^* Ffc(Ai) p fI{X[) e^^* rliX'^)) (3.11) 
k=l 

where the trace is now over the Hilbert space of the component system. 
The last part of the integrand may be written 

TV 

n It (Ffc(A2) e-^* F.iXi) p fI{X[) e^^* fI{X',)) 
k=l 

Tr (Ffc(A2) e-^* Fj,{X,) p ^(A'l) e^^* Fl{X',))f 
gi(Ai+A2-A;-A'2)p(^^ y) + e*(^2-A^)^(^^^ ^) 

_^^i{Xi-K)p(^y^ n) + p{n, n) 

+ei(Ai+A2-A^)^(^^ ^1^^ y) + e^(^2-Al-A^)^(^^ ^|^^ ^) 

+e'^^D{y, n\n, n) + e~^^W{n, n\y, n) (3.12) 

where p(y, y), D{y, y\n, y) etc. denote the probabilities and decoherence functionals 
for each component system, that is 

p{y, y) = IV (P e-^^* PpP e^^*) 



p(n, y) = Tr (^P e-'^* PpP e^^*) 
p{y, n) = Tr (^P e''^* PpP e'^*^ 
p{n, n) = Tr (^P e''^* PpP e^^*) 
D{y, y\n, y)^Tr(p e''^^ PpP e'^^ 



D{n, y\y, y) = TV (p e''^^ PpP e'^^ 

JHt\ 



) 



D{y, n\n, n) = Tr i^P e~^^* PpP e'^^'j 
D{n, n\y, n) = Tv (^P e''^* PpP e'^*^ 
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3.13) 
3.14) 
3.15) 
3.16) 
3.17) 
3.18) 
3.19) 
3.20) 



The probabilities obey the simple relation 



piVi y) + Piji, y) + p{y, n) + p(n, n) = 1 



(3.21) 



The off-diagonal terms of the decoherence functional satisfy 



D{y, y\n, y) + D{y, n|n, n) = 
y\y, y) + D{n, n\y, n) = 



(3.22) 
(3.23) 



The latter relations are particular to decoherence functionals of histories character- 
ized by projections at two moments of time. They are consistent with the general 
property that the sum of off-diagonal terms is zero. Also, since the second relation 
is the complex conjugate of the first, they imply that all the off-diagonal terms of 
the decoherence functional may be recovered from a single complex number, which 
may be taken, for example, to be D{y,y\n,y). 



The integral (3.11), with (3.12) inserted may be evaluate exactly by muliple 
use of the binomial expansion. This is carried out in Appendix A. More useful is 
an approximate evaluation for large N. To do this, note that 



Tr (Ffc(A2) e-'^>^' F,(Ai) fI{\[) e'^^' fI{\'^)) | 

< TV (Ffc(A2) e-'H^' F,(Ai) ^(Ai) e'^^' pliX^)) 
X IV (f,(A'2) e-*^'=* Ffc(A'i) pk Ft(A'i) e^^^' fI{X',)) 



The first relation, the inequality, is an elementary generalization of the inequality 



obeyed by the decoherence functional [10]. The second relation, equality with unity, 
follows from the fact that Fj.{X) obeys the relation 



3(C). Evaluation for Large N 



= 1 



(3.24) 



D(a,a')\'^ < D(a,a) D(a',a') 



(3.25) 



Fk{X)Fl{X) = l 



(3.26) 



as is easily shown. 



16 



The norm of the term raised to the power N when Eq.(3.12) is inserted in 
the integral (3.11) is therefore less than or equal to 1, with equality if and only 
if all the A's are zero. It follows that for very large N, the integral over the A's 
in Eq.(3.11) is dominated by values of the A's close to zero, and the integral may 
therefore be evaluated by exanding about this point. (Note that we could not have 
used the inequality (3.24) if we had made the simplification (Pn2) = Pn2 the 
final projection). 

We now make the following change of variables: 

,^ = ^^, «i = ^. w = ^, 6 = ^. (3.27) 

The integral expression for the decoherence functional then becomes, 

X e^^(^i+^2)p(^^ ^) + e'^^^^p{n, y) + e^^^^p{y, n) + p{n, n) 
+e'(^'^+^^+'^^^^D{y,y\n, y) + e'^-'^^^^^^'^^^^ D{n,y\y, y) 

+e*(/^i+?i)i:)(y, n|n, n) + e'(-^i+^i)i:>(n, n|y, n)] ^ (3.28) 

Since the integrand is independent of /X2) the integral over //2 be carried out. 
As explained above, for very large A^, the integral over the remaining variables 
may then be carried out by expanding the integrand about jJ^i — ii — ^2 = 0- To 
quadratic order, the term in brackets raised to the power N has the form, 

iN 



1 + zx • V - :iJ'M:x. 



(3.29) 



where x is the three- vector (/Ui, ^1,^2)7 v is a three- vector, and M is a symmetric 
3x3 matrix. The components of v may be read off from the following: 

X- V =2 {j){y,y)+p{y,n))ii 

+ 2 {p{y, y) + p(n, y) + D{y, y\n, y) + D{n, y\y, y))^2 
=2po^l + 2pti2 (3.30) 

where the properties of the decoherence functional (3.22), (3.23) have been used, 
and pq and pt denote the probabilities of the alternative "yes" at times and t 
respectively, i.e, 

Pq = Tt{Pp), = Tr (Pe^%e-*^*) (3.31) 
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We also denote the probabilites of the "no" alternatives at times and t by po 
and pt respectively. Hence we have po + po = 1, and pt + Pt = 1- Similarly, the 
components of M may be read off from 

x^Mx =2pod + 2 (My, y) + D{y, y\n, y) + D{n, y\y, y)) ^1^2 

+ 2ptil + 2 {D{y, y\n, y) - D{n, y\y, y)) (3.32) 

Again the properties (3.22), (3.23) have been used. 
For large N, 

(l + ^)^^e^(^-^^') (3.33) 
hence the decoherence functional now has the form 

B(ni,«2l«i,n2) =/rfW</ei<J«2 expOw(*>i - 

-i^l(ni + n[ - 2Npq) - ^6(2^2 - 2Npt) - x^^xj (3.34) 

where here, and in what follows, we will drop overall normalization factors (these 
are readily recovered if required). ^ is a symmetric 3x3 matrix, defined by 

x^Ax = iVx^Mx - -iV(x • (3.35) 

Its explicit components, which will be important below, are 

Aoo = 0, Aoi = (3.36) 

Ao2 = iNIm{D{y,y\n,y)) (3.37) 

All = 2A^P0P0 (3.38) 

Ai2 = 2N [p{y, y) - pm + Re {D{y, y\n, y))] (3.39) 

A22 = 2NptPt (3.40) 

The integrals may then be carried out, with the result, 

/ f / 2 / 2 

-D(ni,n2|ni,n2) = exp \ -a{ni - Ui) - /3(ni + % - 2Np{yo)) 

-i^{ni - n'i){n2 - Npt) 

-w{ni - n'i){ni + n'l - 2Npq) (3.41) 
where the real coefficients a, (3,^,1/ are given by 

a = . .1^ 3.42 

4^1l(z^02)2 
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2^11(^02) 

This is now the decoherence functional for precisely specified values of n. 

To complete the calcuation of the decoherence functional a further coarse- 
graining over n is required, corresponding to imprecise specification of n. This 
involves summing ni, n'^ and n2 over ranges of values, and is most easily achieved, 
at least approximately, by taking n to be continuous and integrating with a Gaus- 
sian smearing. Denoting the coarse-grained values of n by N, we have 



(3.45) 



D{Ni,N2\N[,N2) = J dni dn[ dn2 D{ni,n2\ni,n2) 



.__ exp I J^'-f')^ - - ](3.46) 

(27r(72)3/2 2(72 2(7^ 2(72 ' 

where the ranges of integration are —00 to -\-oo. The coarse-grained variables N 
are understood to have significance only to up order cr. 

The integrals are most easily carried out by changing variables to ni + n'^ and 
^1 ~ ^'i 1 cind by making use of the formula 



/ 



dx 



exp {^a{x — xif' — b{x — ^2)^ + icx^ 



.2 



, ^ ) exp ( -—r- — 7T - ^—77(^1 ~ ^2)^ + i-r — —rrici'Xi + 6x2) 1(3.47) 
Xa + b)J ' \ 4:{a + b) (a + 6) ^ ^ {a + b)^ ^ 

One thus obtains 

D{Ni,N2\N[,N2) 

= exp ( -di{Ni - N[f - /3{Ni + N[- 2Npof 

-e{N2 - Nptf - 0(iV2 - Npt){Ni + N[- 2Npo) 

-z7(iVi - N[){N2 - Npt) - ii>{Ni - N[){Ni + N[- 2Npo)) (3.48) 



where 



^ 2 2 2 2 

^ = T, r"o7T) b^ a + -I- ^ ^ on\ (3.49) 

(1 + 4(726)' 2 (1 + 4(72/?) ^ ^ 
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^3 ^ ^ 

" (1+4^2/}) + (1+4,t26)(1+4,t2/3)2 

~ (l + 4a2^,)(i+4^2^) (3.52) 

= (l + 4^726)(l + 4^72/?) (^-^"^^ 

Eq.(3.48) is the ciccoherence functional for densities of the "yes" ahernative, coarse- 
grained to a width a, and for histories characterized by projections at two moments 
of time. 



3(D). The Degree of Decoherence 



A reasonable measure of approximate decoherence is the size of the off-diagonal 
terms in comparison to the probabilities, Eq.(1.8). Denoting the diagonal parts of 
the decoherence functional hy p{Ni,N2), we have, 

—-_ — _ 1 _ 1 =exp(-(a-/3)(iVi-iVi) I (3.55) 

p(Ari,Ar2)2p(Ar;,Ar2)2 V ) 

hence the degree of approximate decoherence is controled by Q—(i. Since the coarse- 
grained density N has significance only up to order o", the degree of decoherence, 
which we denote e, is at worst given by 

e = exp(|-(a-/3)(j2) (3.56) 

For small cr, o: ~ a and /3 ~ /?. We will see in the next section that, for the 
models we are interested in, the probabilities po,pt, etc. for the component systems 
are of order 1, whilst the off-diagonal terms of the decoherence functional are much 
smaller. ^From Eqs.(3.46)-(3.40), this means that a » the important term in 
a is the term Aq2, and the terms An, A22 and A12 are important only in their 
AT-dependence. It is also convenient to write the coarse-graining parameter cr as a 
fraction / of the total particle number N, so a = fN. The degree of decoherence 
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is therefore of the form 



ivr/2 



(3.57) 



e ~ exp 



{lmD{y,y\n,y)) 



where 



r = PtPt - 



{p{y, y) - Pm + ^^D{y, y\n, y)f 



(3.58) 



This is the expected result, and the main technical result of this paper. The 
degree of decoherence improves with increasing N. It also improves as the degree of 
decoherence of the component systems improves (ie., as \ D{y^y\n,y) \ gets smaller). 
Moreover, decoherence also relies on the factor / not being too small. This means 
that the density n must be partitioned into macroscopically distinct sets for there 
to be sufficient decoherence. 

Note that Eq.(3.57) invites one to use the ratio of (ImL>(y, |/|n, y))^ to F as a 
measure of approximate decoherence of the component systems (rather than, for 
example, the condition (1.8)) - not an immediatey obvious measure. Having said 
that, it will in fact turn out that for the model of the next section, the probabilities 
are of order 1, so F is of order 1. (See figure 3.) 

Note also that the decoherence functional (3.48) involves only the imaginary 
part of the decoherence functional of the component systems. As we shall see in the 
next section, it is possible to have ReZ)(y, y) = but ImD(j/, y\n, y) ^ 0. This 
means that the component systems could be exactly consistent but the total system 
not exactly consistent. There is no contradiction since the decoherence functional 
for the total system is a sum of products of the decoherence functionals for the com- 
ponent systems, so ReZ)(y, y|n, y) = does not imply that ReD(iVi, iV2|-^i) -^2) is 
exactly diagonal. Exact decoherence of the component systems (rather than just 
consistency), however, does imply exact decoherence of the total system. The sig- 
nificance of this, if anything, is to underscore decoherence, Eq.(1.7), as a physically 
more meaningful condition than consistency, Eq.(1.6). 

Given decoherence we may now examine the probabilities for the histories. 
These are given by 




(3.59) 
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This means that the probabihties for Ni only, or for A^2 only, are peaked about 
NpQ or Npt, as one might expect. Beyond this, the form of the probabilitiy distri- 
bution (3.59) is not of much significance. The model is too simple for us to expect 
approximate evolution equations for N. 

There is nothing in this section that refers to the notion of conservation. The 
above result is essentially a combinatoric one. As we will see in the next section, the 
role of conservation is to ensure that \D{y,y\n, y) \ is small, and moreover, becomes 
smaller as the coarse- graining volume increases. 
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4. THE SPIN SYSTEM MODEL 



We now describe a particular model which is of the type discussed in Section 
III, and in which we expect decoherence through approximate conservation. 

4(A). The Model 

The model consists of a chain of very large number M of locally coupled 
"atoms", each of which can be in one of two states, call them spin up and spin 
down. (We follow Ref . [20] ) . The spins interact via the Hamiltonian 

M 

H=-lY.^n- (4.1) 
n=l 

where is the 3- vector whose components are the Pauli matrices. We impose 
periodic boundary conditions, V^v+i = ^l- Up to an additive constant, the 
Hamiltonian may be written 

M 

H = -xYl f''"'"^' (4-2) 
n=l 

Here, p"'^'"'+^ is a spin exchange operator. It has the effect of leaving aligned spins 
alone, and interchanging oppositely aligned spins. For example, 

P^'^ ITi) = lit) (4.3) 

The ground state is the state with all spins pointing up and has eigenvalue 
zero. The next excited states consist of the M states for which one spin is down 
and the rest up, and we denote them 10^), where 

l<^i) = liTT---T) (4.4) 
l02) = ITiT---T) (4.5) 

and so on, i.e., \(f)]^) denotes the state in which spin number k is down and the rest 
are up. These are the so-called spin waves [20]. Then it is readily seen that 

\cf>k) = Id^k+l) (4.6) 
\M = M (4.7) 
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and 

pkMl 10^) = (4.8) 

if n < A; or n > A; + 1. 

The states are not eigenstates of the Hamiltonian, but the combinations 

, , V 1 / 27ri£k\ , , , , , , 

are, with eigenvalues E£ = —2% cos , and £ = 1, • • • M. Their normahzation 

follows from the identity, 

J^f exp(H^!^)=.„„ (4.10) 

fc=l ^ ^ 

where (J^m is the Kronecker delta. 

4(B). Coarse Grainings 

There are two conserved quantities for this model: the total energy and the 
total spin. We will study coarse-grainings which ask for the total spin in a region 
of the chain. We take the crudest coarse-graining, which is to divide the chain into 
two regions, region 1 and region 2, of size Mi and M2, where Mi + M2 = M. We 
will construct projections which ask for the total number of down spins in region 
1. 

The model so far is not yet exactly of the type described in Section III because 
there, a non-interacting Hamiltonian was assumed, whereas here, we have an in- 
teracting one. To get around this, and hence to use the results of Section III, we 
do the following. First, we restrict attention to the subspace in which there is a 
fixed number of spins pointing down. This subspace is invariant under Hamil- 
tonian evolution, so the set of downward pointing spins may be regarded as a 
set of interacting "particles". Second, and more importantly, we assume that 
1 << << M. This is a kind of "dilute gas" assumption, and we expect it to 
allow us to neglect the interaction between downward pointing spins. 

More precisely, we replace the N particle Hilbert space Hjsf with a tensor 
product of N Hilbert spaces Hi, where Hi denotes the Hilbert space of states with 
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one spin pointing down described above. We also take the Hamiltonian Hj^ for the 
N particles to be of the form 

i^jV = ® 1 ® 1 ® • • • 
+ 1 (g) i^l (8) 1 (8) • • • 

+ ••• (4.11) 

where Hi is the Hamiltonian for the states with one spin pointing down, operating 
on Hi- The Hamiltonian is thus of the form required for the results of Section HI 
to be applicable. 

Clearly some work is required to fully justify this approximation. The Hilbert 
space H n does not have exactly the same dimension is the tensor product of the N 
spaces Hi, although they are close for sufficiently large N and M. Also, one would 
expect the effective Hamiltonian to include a term in the Hamiltonian preventing 
two down spins from occupying the same site, and a term describing the interacting 
between neighbouring down spins. What we are assuming, in effect, is that these 
extra terms can be neglected if the gas of downward pointing spins is sufficiently 
dilute. 

Even if this approximation cannot be justified, we could always postulate a 
model of the above type. This would admittedly be unphysical, but at least it 
provides a framework in which we can investigate the mathematical properties we 
are interested in, which is the main aim of this paper. 

The projection Pn onto the total spin in region 1, is constructed from the 
projection for the individual particles, as described in Section HI, Eq.(3.5). For 
each individual particle we introduce the projections. 

Ml 

P=Y1 \'^k)M (4.12) 
k=l 

M 

P= J2 \^k){4>k\ (4.13) 
A:=Mi+l 

In the 1-particle subspace, P and P ask whether the downward spin is in region 1 
or not in region 1 {i.e., in region 2), respectively. We also denote these alternatives 
y and n, respectively, in accordance with the notation of Section HI. As shown in 
Section HI, the computation of the decoherence functional onto histories in which 
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the total spin in region 1 is specified reduces to a computation of the decoherence 
functional for the component systems. 



4(C). The Decoherence Functional of the Component Systems 



The Hamiltonian for the 1-particle subspace TCi may be written, 

M 

Hi = Y,Ei (4-14) 

and hence the unitary evolution operator is 

M 

e.=i 

The Hamiltonian may also be written, 

M 

Hi = -xY.(M{<l^k-l\ + \4>k-l){4>k\) (4.16) 
k=l 

The decoherence functional for a simple two-time history is 

D(ai,a2|«'l,«2) = IV (^Pa,{t)Pa,pPa[) (4-17) 

where Pa2{i) = e^^i^P^ge"*^!*, and Pa denotes P and P. Take an initial state 
which is a superposition states in which the particles is in region 1 and 2, 

I*) = ^ {\4>k,) + (4-18) 
where 1 < A;i < Mi, and Mi + 1 <k2<M. Then 

Pa\^)mPa' = l\cl>kJ{<l>kJ (4.19) 

where for the moment we let a take the values 1,2, corresponding to y,n. The 
decoherence functional then is 

D{ai,a2\a[,a2) = ^{(t>k^jPa2{t)\(f>k„,) (4.20) 
For this, we then have, for example 

D{y,y\n,y)=^{<Pk^\P{t)\<Pk,) (4.21) 

P{n,y) = l{4>k,\Pmk,) (4.22) 
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and so on. 

Let = P- Then it is readily shown that 

M M 



where 

,(,,)^E-P(-?^^) (4.24) 

k=l ^ ^ 

Now, writing x = 27r{£ — £')/M, we have 

^ y e-^^^ = — ^ (4.25) 

The properties of the decoherence functional may be understood through d{£,£'). 
Clearly, if Mi = M, then d{£,£') = MSn,, and P{t) would be diagonal in the |0n)'s. 
More generally, if Mi and M are very large, then d{£, £') is very small when £ ^ 
and d{£,£) = Mi. 

Curiously, when Mi = M2 the real part of the decoherence functional vanishes. 
Thus, this coarse-graining exhibits exact consistency, but not exact decoherence. 
It seems reasonable that this sort of consistency, where the imaginary part of the 
decoherence functional is non-zero, will only occur due to symmetries in the system 
and choice of histories. 



4(D). Numerical Results 



The decoherence functional for the component systems may be computed nu- 
merically using Eq.(4.20). There are three parameters, t, M and Mi. The deco- 
herence functional for various ranges of the parameters is computed below and the 
results plotted in figures 1-3. 



4(E). Summary of Results 



The results of the numerical calculation may be concisely summarized as fol- 
lows. The probabilities j/), p{y,n), etc. are of order 1. The off'-diagonal terms 
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of the decoherence functional, divided by the probabihties, are typically of or- 
der M~2, hence the degree of decoherence for the 1-particle case is of this order. 
This degree of decoherence is not particularly good (compared, for example, to the 
quantum Brownian motion models, where one typically finds that the degree of 
decoherence is an exponential function of the coarse-graining parameters [10]). It 
is no better than one would expect from the overlap of two arbitrary vectors in a 
Hilbert space of dimension M [21]. Physically, this is not surprising, because the 
histories in question differ by just one quantum of spin, and one would not expect 
the interference between such spins to be substantially supressed. 

However, when inserted into the expression for the decoherence functional of 
the A^-particle case, (3.48), a vastly improved degree of decoherence is obtained, for 
sufficiently large A^. Moreover, the important feature about the decoherence func- 
tional (4.20) is that the off-diagonal terms decrease as the size Mi of the "smearing 
volume" increases from zero (until Mi reaches the value about M/2). Correspond- 
ingly, the off-diagonal terms of the decoherence functional of the A'"-particle case 
(3.48) decrease very rapidly as Mi increases from zero. 



28 



5. SUMMARY AND CONCLUSIONS 



This paper marks a first step towards the problem of computing the decoherence 
functional for hydrodynamic histories. We considered systems of essentially non- 
interacting particles and studied histories of projections onto densities of those 
particles. Our principal aim was to show, in the context of some simple models, 
how decoherence can come about as a result of approximate conservation, in tune 
with the general ideas put forward by Gell-Mann and Hartle [4,5]. 

The results of Section IV show explicitly how approximate decoherence is re- 
lated to approximate conservation: the degree of decoherence increases as the 
smearing volume increases (at least, for volumes less than half the total volume of 
the system). The degree of decoherence, however, is not very good. This is because 
we considered single particle models in Section IV, and the histories differ by just 
one quantum of spin. The histories are therefore not "macroscopically distinct" 
and one would expect that interference between them could still be quite notice- 
able. Differently put, it is because the fluctuations in variables in question {e.g., 
number density) are comparable to the values of the variables themselves. 

To obtain adequate decoherence, it is necessary to have a large number N 
of particles and partition the number density by large ranges a. It is here that 
the results of Section III came in. We found that when the histories do differ 
by macroscopically significant amounts, interference is destroyed very efficiently. 
Clearly what is happening here is that the fluctuations in the variables are much 
smaller than the variables themselves. 

^From this we conclude that, in these simple models, decoherence requires two 
distinct phenomena: approximate conservation, and large particle number parti- 
tioned into large ranges. Approximate conservation ensures that the dynamical 
variables projected onto become sufficiently slowly varying for sufficiently large 
smearing volume. Large particle number partitioned into large ranges ensures that 
the quantum fluctuations in the local densities are smaller than the values of the 
variables themselves. These conclusions concure with the general expectations ex- 
pressed in Refs.[4,5]. 

The restriction to non-interacting particles may seem rather unrealistic. It 
means that conservation may already be seen at the one particle level, whilst for 
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an interacting theory, it is only seen for the whole collection of particles. Although 
physically unrealistic, it had the mathematical advantage that the effects of large 
particle number and approximate conservation could be cleanly seperated. More- 
over, one would not expect the inclusion of interactions to substantially modify 
our conclusions (although this is clearly an important extension to carry out). The 
point is that here, it is the approximate conservation and large particle number 
that produce decoherence, not the interactions. This is in stark contrast to the 
vast majority of other models studied in the literature on decoherence, in which it 
is interactions (usually with another system) that are held responsible for decoher- 
ence. 

This paper is, as stated, a first step in an extensive investigation, and there 
are therefore many ways in which it may be developed. Perhaps the next step is 
to seek a more sophisticated model in which one would expect the probabilities 
for histories to be peaked about interesting approximately deterministic evolution 
equations. 

A technically very different but similar in spirit investigation is that of Calzetta 
and Hu [11], who considered the decoherence of histories characterized by the values 
of n-point functions for fields. They did not make any contact with the notion of 
approximate conservation, but it is similar to our work in that it is, to the best of 
our knowledge, the only other concrete calculation of decoherence which avoids a 
system-environment split. It would be of interest to find more detailed connections 
between their work and ours. 

These and other questions will be the topic of future publications. 
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APPENDIX A: EXACT EVALUATION 
OF THE DECOHERENCE FUNCTIONAL 



The integrals over the A's may be evaluated exactly, and we present here the 
result of that calculation. 



= dXi dX[ j"" d\2 e-2^^i^i+2iA;n;-2iA2n2 

_^g2i(Ai-A;)^^^^ n) + pin, n) 

^^HM+X2)D{y^ y|n, y) + e2^(^2-^'i)l)(n, y\y, y) 
+e^^^'^D{y, n\n, n) + e~'^^^W{n, n\y, n) 

Now note that the term in brackets raised to the power A*" may be written 

N 



A + Be 



2i\2 



k=0 ^ 



N-k j^k g2iA2fc 



using the binomial expansion, where 

A = e^*(^i~-^'i n) + p{n, n) + e^^^W{y, n\n, n) + e~'^^^W{n, n\y, n) 
B = e2^(^i-^i)p(y, y) + p{n, y) + e^'^^D{y, y|n, y) + e-'^'^W{n, y\y, y) 

The integral over A2 may now be evaluated, with the result, 

{2TTy J \n2j 
We now use the binomial expansion again on ^^""'2 ^^ite 

A = e-2*^ic + D 



where 



B = e'^'^'^E + F 



C = e^*^ip(y, n) + D{n, n\y, n) 
D — pin, n) + e^^'^^Diy, n\n, n) 
E = e'^'^^p{y,y) + D{n,y\y,y) 
F = p{n,y) + e^'^W{y,y\n,y) 
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(Al) 



(A2) 
{A3) 

{A5) 



(AG) 

(AS) 

(AlO) 
(All) 



Expanding A-^ B'"''^, the integral over X'^ may then be carried out, with the 

result, 

N—n2 m 



k+j,n[ 



iAl2) 



Repeating the binomial expansion of etc. on final time, the last integration 
over Ai may be performed, and we arrive at the result, 

N-n2 n2 k N-fi2-k j fi2-j 



D{ni,n2\n'i,n2) 



N 

m 



X 



X 



X 



X 



X 



E EE E EE 

fc=0 j=0 £=0 m=0 r=0 s=0 

'A^ — f n2 



k 



^k+j,n[ ^i+m+r+s,ni 



N -n2-k 
m 

r 



[p(n,n)]^-^2-^-'« [D{y,n\n,n)r 
[D{n,y\y,y)y~'' \p{y,y)Y 

\p{n,y)f^-^-' [D{y,y\n,y)Y (A13) 



For reasonably small A^, this expression can be evaluated numerically. For large 
(and ni, n'^, ^2), the sums can be estimated by replacing them by integrals, and 
assuming that the dominant contributions come from the middle of the ranges, i.e., 
from values of j, k,£,m,r,s that are large (although we have not carried this out 
explicitly). 
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Figure 1. The probabilities p{l, 1), p{l, 2), p{2, 1), p{2, 2) as a function of the coarse- 
graining Ml. This plot was produced at t = 1000 for a string of M = 1000 spins; the 
energy constant is x = 1- Note that for Mi close to 500 (half the spins), the probabilities 
are all close to 0.25. 

Figure 2. \D[1, 1;2, 1] |/ a/p(1, 1)p(2, 1) vs. Mi. For this plot the parameters are t = 
1000, M = 1000, X = 1. Note this quantitity is highly peaked for low values of Mi, 
indicating that better decoherence results from coarser graining. 

Figure 3. |Imi:)[l, 1; 2, 1] p/F vs. Mi. For this plot the parameters are t = 1000, 
M = 1000, X = 1- This also exhibits peaking for low values of Mi, indicating good 
decoherence for coarse graining. Note also that the absolute value of this quantity is low; 
this will produce very good decoherence for many-spin systems, as shown in equation 
(3.57). 



36 



5e-05 



4.5e-05 - 




3e-05 - 

2.5e-05 - 

2e-05 - 

1.5e-05 - 

1e-05 - 

5e-06 - 








IImD(l,l;2,l) |2/r vs. Ml, t= 1000 
1 1 1 



100 200 300 400 



500 
M1 




600 700 800 900 1000 



